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In this article we consider two spin—1/2 chains described, respectively, by the thermodynamic 
limit of the XY model with the usual two site interaction, and an extension of this model (without 
taking the thermodynamics limit), called XYT, were a three site interaction term is presented. 
To investigate the critical behaviour of such systems we employ tools from quantum information 
theory. Specifically, we show that the local quantum uncertainty, a quantity introduced in order 
to quantify the minimum quantum share of the variance of a local measurement, can be used to 
indicate quantum phase transitions presented by these models at zero temperature. Due to the 
connection of this quantity with the quantum Fisher information, the results presented here may be 
relevant for quantum metrology and quantum thermodynamics. 


I. INTRODUCTION 

Characterized by a change in the macroscopic order 
of a physical system, phase transitions are ubiquitous in 
nature. Contrary to thermodynamical phase transitions, 
which occur at finite temperature, quantum phase tran¬ 
sitions (QPTs) take place at zero temperature, driven by 
quantum fluctuations [T]. Indeed, being linked with ener¬ 
gies level crossings in the ground state, QPTs can occur 
even at finite temperature, as long as the system can¬ 
not be excited by the thermal fluctuations. A candidate 
system to study these kinds of phenomena is a quantum 
spin chain, which can present several distinct phases (Tj. 

Many interesting physical systems exhibit QPTs, 
among which we mention the superfluid to insulating [2] 
and the superconducting to insulating [3] transitions for 
bosonic systems while, for fermionic ones, we have the 
spin wave density to a magnetic phase 0 Ej and the 
quantum spin glass to paramagnetic phase transitions [5] 
as examples. A remarkable feature of QPTs is that, in 
the neighbourhood of a critical point, the system becomes 
strongly correlated, and its correlation length diverges at 
this point 0 . This has suggested that quantum correla¬ 
tions should be maximum near the transition point, fact 
that was verified in seminal works about entanglement 
measures in critical systems (see Ref. [5] for a re¬ 
view regarding entanglement in many-body systems). As 
a natural consequence, other measures of quantum corre¬ 
lations were also employed to investigate the properties 
of QPTs, such as quantum discord [MIS] and related 
quantities [MM 

Within the last decade, quantum correlations other 
than entanglement were discovered [20] El] and their im¬ 
portance for quantum information processing has only 
increased [2121 - E5] . The origin of such correlations lies in 
the existence of quantum coherence among different par¬ 
titions of a quantum system, being present even in sep¬ 
arable states (see Ref. |25] for a recent review). Among 
all the quantum correlation measures, the so-called local 
quantum uncertainty (LQU) stands out as a bona fide 
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measure with a good computability and a close mathe¬ 
matical expression for any qubit-qudit bipartite system 
EZj. This measure is based on the skew information, 
introduced in Ref. [2&j to quantify the uncertainty in 
the measurement of an observable. The LQU is defined 
just as the minimazation of this quantity over all pos¬ 
sible observables. Beyond its importance as a quantum 
correlation quantifier, LQU is also related with the quan¬ 
tum Fisher information )29l4-> 1 1 . which makes it very im¬ 
portant for quantum metrology. Moreover, it has been 
shown that it is also connected with the speed of evolu¬ 
tion of a quantum system m These are some of the 
reasons that lead us to choose LQU as a tool to investi¬ 
gate critical behaviour. 

Here we study local quantum uncertainty in critical 
systems considering two spin-1/2 chains, described by 
the anisotropic XY Hamiltonian and the XY model 
with three spins interaction, named XYT model [M 
134] . which contains the spin-liquid phase [35; that only 
recently could be experimentally investigated [3B]. We 
consider the zero temperature case, making some com¬ 
ments about finite temperature in the last section. 

Several works have addressed the critical behaviour 
of the XY model from the point of view of quantum 
information, considering quantum and classical correla¬ 
tions quantifiers as tools for the identification of quantum 
phase transitions (see, e.g., Refs. nm i2si and references 
therein). Recently, the local quantum uncertainty was 
also computed for some instances of this model and its 
usefulness for the identification of the QPT was verified 
m■ However, the case of the XYT spin chain received 
much little attention. In Ref. [34] the authors demon¬ 
strated that quantum discord [20] and classical correla¬ 
tions [3Ij are good indicators for the critical behaviour of 
this model, while entanglement measures were considered 
in Ref. 2S] ■ The purpose of the present article is to con¬ 
tribute to this debate by investigating how local quantum 
uncertainty behaves for the XYT model. We also provide 
a complementary analysis of the usual XY Hamiltonian, 
presented in Ref. m- While the XY Hamiltonian is 
investigated in the thermodynamic system, for the XYT 
model we consider a finite chain. 

The paper is structured as follows. In Section II we 
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review the definition of the local quantum uncertainty 
while Section III is devoted to present both models and 
our results. The final considerations are done in IV. 


II. LOCAL QUANTUM UNCERTAINTY 


state cannot be an eigenstate of any local observable, thus 
inevitably leading to a non-vanishing variance. 

In the following section we employ this tool to study 
the critical behaviour of two spin chains, showing that 
it is a very good candidate for pinpoint quantum phase 
transitions in both models. 


The LQU was introduced in order to quantify the min¬ 
imal quantum uncertainty that could be attainable in a 
measurement of a local observable E3- This minimal 
uncertainty comes from the existence of quantum corre¬ 
lations between the system being measured and the rest 
of the universe. Based on this fact, a measure of quan¬ 
tum correlations was then proposed m- Let p be the 
state of a bipartite system and I\ some local observable 
on one of the partitions, from here on called the system. 
LQU is defined as 

U(p) = minl(p, / ® A'), (1) 


where 

T(p,/®A) = ~Tr ([^,A] 2 ) 

is the so-called skew information [581EZ2 1 which is a mea¬ 
sure of the non-commutativity between the state and the 
observable. The minimum is taken over the set of all 
observables acting on the considered partition. Strictly 
speaking, the definition of IA should take into account 
the spectrum of the observable being measured. Each 
different choice of the spectrum leads us to a different 
family of measures. However, if the system is a qubit, 
as in the present article, the spectrum does not change 
the behaviour of U (except for a global multiplicative 
constant), therefore we do not explicitly considered it in 
the definition (jT]) , as originally done m- Moreover, as 
we are dealing only with qubits (spin— 1/2 particles), we 
consider the following general form for the local observ¬ 
able K = r • <x, where a is the Pauli vector operator and 
|r| = 1. This enables one to compute a closed, analytical, 
form for LQU j27( 

u{p) = 1 - X max {W}, 

where the matrix W has the following elements 

( W)ij = TV {s/p(ai ® 1 l)s/p((Jj <g> I)}, 

and Xmax denotes the maximum eigenvalue of the matrix 
W. 

As LQU quantifies the minimum quantum uncertainty 
in a measurement performed on the system, it is zero only 
when the state of the system commutes with some observ¬ 
able, which implies that the local system is not quantum 
correlated with anything else. Note that the existence of 
this uncertainty (even in the presence of a flawless appa¬ 
ratus) is caused only by the probabilistic character of the 
quantum measurement. In other words, if the system of 
interest is correlated with some other system, its reduced 


III. MODEL AND RESULTS 

As stated in the Introduction, we consider two distinct 
models, the AT spin chain with two and three spins in¬ 
teraction. The case of two spins interaction was recently 
considered in Ref. m■ However, the focus of that work 
is the so called local quantum coherence, as measured by 
the skew information defined in the preceding section, in 
the XY model. Although U is an extremization of X, 
here we extend the analysis reported in Ref. [37] con¬ 
sidering the entire range of the anisotropic parameter. 
Moreover, we focus in the case of three spins interaction 
model, which was not considered in Ref. m- 


A. The XY model with two spins interaction 


We describe here the paradigmatic example of the 
XY model, which govern the dynamics of a spin—1/2 
chain under the effect of a magnetic field in z—direction 
anisotropically interacting in the xy plane. The Hamil¬ 
tonian of this model is given by 

1 N 

H = ~ 2 51 [C 1 +' K a 'i ■ i + (! - + Xa iP) 

3 =1 

with 7 £ [0,1] and N being, respectively, the anisotropy 
parameter and the dimension of the chain, and A stands 
for the strength of the external field. For 7 = 0 the 
Hamiltonian reduces to the XX model while for 7 = 1 
we obtain the Ising Hamiltonian in a transverse field. For 
any other value (0 < 7 < 1) H belongs to the Ising uni¬ 
versality class. A is the reciprocal of the external trans¬ 
verse magnetic field and crj' is the i —th Pauli matrix for 
spin j. 

Through Jordan-Wigner and Bogoliubov transforma¬ 
tions, the Hamiltonian ([ 2 ]) can be diagonalized, in the 
thermodynamic limit (N —> 00 ), in terms of spinless 
fermions operators [T]. Remembering that the model 
posses translational and parity symmetries [8j, the re¬ 
duced density matrix of two arbitrary spins can be writ¬ 
ten as 



Ion + W z ){v'q + &n) 


k =1 


'k_k\_k_k 

\ a 0 a n) a 0 a n 


(3) 


The details of the correlation functions appearing in the 
above equation are given in Appendix A. It is important 
to note here that we are dealing with the so-called ther¬ 
mal ground state, which is the zero temperature limit of 
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the canonical ensemble. Moreover, we are disregarding 
the effects of symmetry breaking (see Ref. [13] to further 
details regarding this issue). 

This model presents a second order quantum phase 
transition at the critical point A c = 1, for all values of 7 , 
which is of an order-disorder type, separating the ferro¬ 
magnetic from the paramagnetic phases [T]. 

Figure [l] shows the local quantum uncertainty and its 
derivative for three relevant cases for the XY model, con¬ 
sidering first-neighbours interactions. As we can see, the 
QPT at the critical point A = 1 can be clearly identified 
by IA , as confirmed by its first derivative. The more pro¬ 
nounced identification is obtained by the XX model. We 
also observe that, beyond this divergence in the deriva¬ 
tive of IA at the critical point, other discontinuities points 
are present. These points are related with the minimiza¬ 
tion process involved in the definition of U and not with 
a QPT, which is related with a divergence presented in 
the derivative of the density matrix elements, as occurs 
at the critical point. The results presented in this sec¬ 
tion are in complete agreement with the study recently 
reported in Ref. where the local quantum coherence 
(as measured by the Wigner-Yanasse information) and 
the local quantum uncertainty were studied for the XY 
model considering the cases 7 = 0.5 and 7 = 1. 




A 


FIG. 1. Local quantum uncertainty for three distinct in¬ 
stances of the XY model: The XX model (7 = 0), one ex¬ 
ample of the Ising universality class (7 = 0.5) and the Ising 
Hamiltonian in a transverse field (7 = 1). The critical point 
A = 1 can be clearly identified by U, as confirmed by its first 
derivative, specially for the case of the XX model. 


U 
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FIG. 2. Local quantum uncertainty along with its derivative 
as function of A and 7 for the XY model. Although present¬ 
ing different intensities, the critical line A = 1 can be clearly 
identified by the derivative of U for all values of 7 . The in¬ 
tensities of the quantities increase from the blue color to the 
red one. 

To finish our analysis of this model, in Fig. [2] we show 
the behaviour of IA 1 along with its first derivative with 
respect to A, as function of A and 7 . It is clear that the 
critical line A = 1 can be identified by IA in the entire 
range of 7 , as confirmed by its derivative. 

Although the analysis presented in this section were 
carried out considering only first-neighbours interactions, 
we have also performed the same calculations consider¬ 
ing interactions up to the 20 th-neighbour interaction and 
the qualitative results obtained were the same, thus con¬ 
firming the ability of IA to signal the critical behaviour of 
this model. 


B. The XY with three spins interaction 

As stated in the introduction, the XYT model is an ex¬ 
tended XY model with includes a three-spin interaction, 
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governed by the Hamiltonian 

1 N 

H =-2^ K 1 + T'K 1 a J+i + C 1 - ^jOj+i 

3 =1 

+ Ac 7] + a K + aj_ l( rj +1 ) of] , (4) 

with AT, 7 , and A defined in the preceding subsection, 
while a is the intensity of the three spin interaction. In 
the limit a-> 0 we obtain the usual two spin interaction 
XY model. 

Through Jordan-Wigner and Bogoliubov transforma¬ 
tions, the Hamiltonian Q can also be diagonalized in 
terms of spinless fermions operators pLj 

M 

H= J2 2£ fc (4<4-l/2), (5) 

k——M 

where d k satisfies the anti-commutator relations for 
fermions and M = N/2 for even TV or M = ( N — l)/2 
for odd N. e k = \/Cfc + ('ismx k ) 2 , C fc = A - cosx k - 
2a cos 2a : k and x k = 2nk/N. Considering the case of 
global thermal equilibrium state with inverse tempera¬ 
ture j3 (Boltzmann constant equal to one), the reduced 
density matrix of two arbitrary spins in the chain is given 
by Eq. <®> , but with different formulas for the correlation 
functions 1341 . Such formulas are given in Appendix B. 

The critical properties of the XYT model were studied 
in Ref. [32] and are depicted in Fig. [3] where the ferro¬ 
magnetic and two spin-liquid phases can be identified. It 
is interesting to note that these phases are independent 
of 7 , i.e., the phase diagram of Fig. [3] is the same in¬ 
dependently of the value of 7 . Therefore, to study the 
behaviour of the local quantum uncertainty near the crit¬ 
ical points of this model, it is enough to pick just a single 
value of 7 . Of course, the value of U (and of its derivative 
with respect to A) will depend on 7 , but the qualitative 
behaviour, the one which we are interested here, will be 
the same. In other words, if U is able to signal the phase 
transition for a single value of A, it will be for all other 
values, just the intensities will change. 

Figure [4] shows the variation of IT along with its deriva¬ 
tive for the XYT model. To make these figures we chose 
7 = 0.5, N = 2000 and considered the zero temper¬ 
ature case with first-neighbour interaction. As we can 
clearly see from the figure, the first derivative of the local 
quantum uncertainty is able to signal both critical lines 
that appear in Fig. [3| Although not shown here, we 
performed several numerical studies considering different 
chain sizes as well as neighbour interactions other than 
the first. The results show that qualitative behaviour 
depicted in Fig. [4] do not change, only the intensities 
vary, i.e., IT is perfectly able to signal the critical lines 
appearing in Fig. [3] in all cases. 



FIG. 3. Phase diagram for the XYT model where the Ferro¬ 
magnetic and the two spin-liquid phases are shown [32] . 


IV. DISCUSSIONS 

In this paper we had addressed the local quantum un¬ 
certainty in critical systems described by the XY model 
with two and three spins interaction, the so-called XYT 
model. Our results show that the LQU is a very useful 
tool in the investigation of critical behaviour. Through 
divergences of the derivative of U , the critical lines of 
both models could be identified. 

All the analysis reported here were performed consid¬ 
ering the zero temperature case and only first-neighbours 
interactions. However, we also performed extensive nu¬ 
merical calculations for different neighbours, obtaining 
the same qualitative behaviour. The only change was 
in the intensities of the considered quantities, IT and its 
derivative. While the XY model was investigated in the 
thermodynamic limit, for the XYT we considered a fi¬ 
nite spin chain. Again, the size of the chain only changes 
the quantitative character of the results, with the quali¬ 
tative behaviour, i.e., the ability of IT to signal the phase 
transitions of the model, being the same. 

It is important to observe here that it is not always true 
that a sudden change in the behaviour of IT is connected 
with a quantum phase transition. As a consequence of 
the extremization process involved in the definition of 
IT, the optimal observer can change while a Hamiltonian 
parameter is continuously modified, as also occurs with 
other measures of correlations [39] [40] (see also the dis¬ 
cussion in Ref. Ezj). We have also analysed the case of 
finite temperature (see [111 |14j for the quantum discord 
case). However, in this case, the divergence of the deriva¬ 
tive (with respect to A) of the local quantum uncertainty 
appears only as a local maximum at the critical point A c . 
We can then look at these points in order to signal the 
quantum phase transition. 

Due to the connection of IT with the quantum Fisher in¬ 
formation, the results presented here indicates that crit¬ 
ical systems could be explored in the field of quantum 
metrology mi 01. Moreover, the connection with ther- 
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site z —magnetization is given by [ 8 : 



‘> = - 


(1 + A cos </>) tanh(/3 ui^) 
2 TTUlA 


#» 


with Wfe = ( 7 A sin </>) 2 + (1 + A cos </) 2 /2, f3 = 1 /kgT 

and ks is Boltzmann’s constant. The two-site correla¬ 
tions functions are given by 


(«> = 


/ „X\ _ 

\°b <T n) ~ 


Gx 

Go •• 

’ G _„+2 

g 2 

Gx •• 

’ G_„+3 

G n 

G n -1 ' ’ 

• Gx 

G- 1 

G—2 • 

■■ G- n 

G 0 

G-x • 

’ ’ G-n+l 

3 . . 

1 

to 

G n - 3 • 

■■ G-x 


and 


(a 0 V~) = (a 2 ) 2 - G n G- 

where 


G n = f [ C0S ( n ^)(l + A cos <f>)~ 

J 0 27TW0 

7 A sin(n</) sin (/]<T</>. 


(Al) 


FIG. 4. Local quantum uncertainty along with its derivative 
as function of A and a for the XYT model. We consider zero 
temperature, only the first-neighbour interaction, N = 2000 
and 7 = 0.5. 

modynamics [43] can shed fresh light on questions about 
the irreversible entropy production |44] or the work dis¬ 
tribution |i 45] in critical systems. 
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Appendix A: Correlations for the XY model 

As stated in the main text, we consider here the XY 
model in the thermodynamic limit (TV —> 00 ). The one- 


The zero temperature limit, case considered in this paper, 
of these expressions is easily obtained by taking /3 —>• 0. 

Appendix B: Correlation functions for the XYT 
model 

The 2 —magnetization is given by |34j 

(<r z ) = tanh ’ 

k £k 

where TV is the size of the chain and 

£k = \J Cfc +7 2 sin 2 (x k ), 

Ck = A — cos (xfc) — 2a cos (2x fc ) and X/. = 2nk/N. The 
two-site correlation functions are defined as 
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and 


WO = w z ) 2 -g n g- r 


with 


9n = - 


E 


k 


1 

Ne k 


tanh (/?£fc). 


[cos(a;fen)Cfe + 7 sin(a’ fc n) sin(xfc)] x 


Note that the above expressions, for the limit a —► 
0 give us the relevant correlation functions for the XY 
model without taking the thermodynamics limit. As in 
the preceding case, the zero temperature limit is easily 
obtained here. 


[1] S. Sachdev, Quantum phase transitions (Cambridge Uni¬ 
versity Press, Cambridge, 2011). 

[2] M. P. A. Fisher, P. B. Weichman, G. Grinstein and D. S. 
Fisher, Phys. Rev. B, 40 , 546 (1989). 

[3] Y. Liu and A. M. Goldman, Physica D 83, 163 (1995). 

[4] N. D. Mathur, F. M. Grosche, S. R. Julian, I. R. Walker, 
D. M. Freye, R. K. W. Haselwimmer and G. G. Lon- 
zarich, Nature 394, 39 (1998). 

[5] H. V. Lohneysen, A. Rosh, M. Vojta and Wolfle, Rev. 
Mod. Phys. 79 , 1015 (2007). 

[6] S. Sachdev, Phys. World 10 , 25 (1994). 

[7] A. Osterloh, L. Amico, G. Falci and R. Fazio, Nature 
416 , 608 (2002). 

[8] T. J. Osborne and M. A. Nielsen, Phys. Rev. A 66, 
032110 (2002). 

[9] L. Amico, R. Fazio, A. Osterloh and V. Vedral, Rev. 
Mod. Phys. 80 , 517 (2008). 

[10] M. S. Sarandy, Phys. Rev. A 80 , 022108 (2009). 

[11] J. Maziero, H. C. Guzman, L. C. Celeri, M. S. Sarandy 
and R. M. Serra, Phys. Rev. A 82, 012106 (2010). 

[12] J. Maziero, L. C. Celeri, R. M. Serra and M. Sarandy, 
Phys. Lett. A 376 , 1540 (2012). 

[13] M. S. Sarandy, T. R. de Oliveira and L. Amico, Int. J. 
Mod. Phys. B 27 , (2013). 

[14] T. Werlang, C. Trippe, G. A. P. Ribeiro and G. Rigolin, 
Phys. Rev. Lett. 105 , 095702 (2010). 

[15] T. Werlang, G. A. P. Ribeiro and G. Rigolin, Phys. Rev. 
A 83 , 062334 (2011). 

[16] C. C. Rulli and M. S. Sarandy, Phys. Rev. A 84, 042109 

( 2011 ). 

[17] L. Justino and T. R. de Oliveira, Phys. Rev. A 85, 052128 

( 2012 ). 

[18] W. W. Cheng, C. J. Shan, Y. B. Sheng, L. Y. Gong, S. 
M. Zhao and B. Y. Zheng, Physica E 44, 1320 (2012). 

[19] C.-J. Shan, W.-W. Cheng, J.-B. Liu, Y.-S. Cheng and 
T.-K. Liu, Sci. Rep. 4, 4473 (2014). 

[20] H. Ollivier and W. H. Zurek, Phys. Rev. Lett. 88, 017901 

( 2001 ). 

[21] L. Henderson and V. Vedral, J. Phys. A: Math. Gen. 34, 
6899 (2001). 

[22] M. Gu, H. M. Chrzanowski, S. M. Assad, T. Symul, K. 
Modi, T. C. Ralph, V. Vedral and P. K. Lam, Nature 
Phys. 8, 671 (2012) 

[23] B. Dakic, Y. O. Lipp, X. Ma, M. Ringbauer, S. 


Kropastschek, S. Barz, T. Paterek, V. Vedral, A. 
Zeilinger, C. Brukner and P. Walther, Nature Physics 
8 , 666 ( 2012 ). 

[24] A. Datta, A. Shaji and C. M. Caves, Phys. Rev. Lett. 
100 , 050502 (2008). 

[25] B. P. Lanyon, M. Barbieri, M. P. Almeida and A. G. 
White, Phys. Rev. Lett. 101 , 200501 (2008). 

[26] K. Modi, A. Brodutch, H. Cable, T. Paterek and V. Ve¬ 
dral, Rev. Mod. Phys. 84, 1655 (2012). 

[27] D. Girolami, T. Tufarelli, and G. Adesso, Phys. Rev. 
Lett. 110 , 240402 (2013). 

[28] E. P. Wigner and M. M. Yanasse, Proc. Natl. Acad. Sci. 
USA 49, 910 (1963). 

[29] S. Luo, Phys. Rev. Lett. 91, 180403 (2003). 

[30] D. Petz and C. Ghinea, arXiv:1008.2417 (2010). 

[31] S. Luo, Proc. Amer. Math. Soc. 132 , 885 (2003). 

[32] M. Suzuki, Prog. Theor. Phys. 46, 1337 (1971). 

[33] I. Titvinidze and G. I. Japaridze, Eur. Phys. J. B 32 , 383 
(2003). 

[34] Y-C. Li and H-Q. Lin, Phys. Rev. A 83, 052323 (2011). 

[35] P. W. Anderson, Mater. Res. Bull. 8 , 153 (1973). 

[36] L. Balents, Nature 464, 199 (2010). 

[37] G. Karpat, B. Qakmak and F. F. Fanchini, Phys. Rev. B 
90 , 104431 (2014). 

[38] A. A. Zvyagin, Phys. Rev. B 80, 014414 (2009). 

[39] J. Maziero, L. C. Celeri, R. M. Serra and V. Vedral, Phys. 
Rev. A 80 , 044102 (2009). 

[40] R. Auccaise, L. C. Celeri, D. O. Soares-Pinto, E. R. de 
Azevedo, J. Maziero, A. M. Souza, T. J. Bonagamba, R. 
S. Sarthour, I. S. Oliveira and R. M. Serra, Phys. Rev. 
Lett. 107 , 140403 (2011). 

[41] V. Giovannetti, S. Lloyd and L. Maccone, Nature Pho¬ 
tonics 5, 222 (2011). 

]42] U. Marzolino and T. Proseii, Phys. Rev. A 90 , 062130 
(2014). 

[43] K. Micadei, R. M. Serra and L. C. Celeri, Phys. Rev. E 
88 , 062123 (2013). 

[44] S. Deffner and E. Lutz, Phys. Rev. Lett. 107 , 140404 

( 2011 ). 

[45] E. Mascarenhas, H. Braganca, R. Dorner, M. Franca San¬ 
tos, V. Vedral, K. Modi and J. Goold, Phys. Rev. E 89, 
062103 (2014). 



